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Abstract 

This paper deals with the derivation of Eulerian models of cell populations out of a mi- 
croscopic Lagrangian description of the underlying physical particle system. By looking at 
the spatial distribution of cells in terms of time-evolving probability measures, rather than at 
individual cell paths, an ensemble representation of the cell colony is obtained, which can be 
either discrete or continuous according to the spatial structure of the probability. Remark- 
ably, such an approach does not call for any assumption of continuity of the matter, thus 
providing consistency to the same modeling framework across all levels of representation. In 
addition, it is suitable to cope with the often ambiguous translation of microscopic arguments 
into continuous descriptions. Finally, by grounding cell dynamics on cell-cell interactions, it 
enables the concept of multiscale dynamics to be introduced and linked to the sensing ability 
of the cells. 

Keywords: cell populations, probability measures, discrete vs. continuous descriptions, 
multiscale dynamics 
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1 Introduction 

Models describing the spatiotemporal patterns of cell populations play a prominent role in Math- 
ematical Biology, as they underlie plenty of different applications: from the formation of vascular 
plcxcs (angiogenesis, vasculogenesis) to the growth of cell layers and tumor masses, possibly in 
interaction with surrounding tissues (cancer invasion). The level of description ranges usually from 
the microscopic scale of single cells to the macroscopic scale of large cell aggregates and tissues. 

Looking at relatively recent contributions (see e.g., [3J [3J [T3J [T7] and references therein), 
it appears that more and more importance is being given to nonlocal aspects of cell dynamics. 
Experimental investigations have indeed demonstrated that a relevant component in cell behavior 
is the ability of single cells to sense the surroundings and take decisions accordingly, in a non 
strictly mechanical manner. Apart from phenomenological arguments, at the macroscopic level 
nonlocal terms are often hard to justify within the classical paradigms of the continuum approach. 
In fact they generally resort to concepts, such as the sensing neighborhood of cells with related 
sensing radius, which actually pertain to the length scale of a single cell and are therefore not 
fully compatible with a macroscopic scale viewpoint. Microscopic models are conceptually more 
appropriate, however the problem still remains of how to translate microscopic ideas into a repre- 
sentation of the system not necessarily focused on individual cells. Indeed, continuous descriptions 
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are often useful over discrete agent-based ones, especially in case of large cell colonies, for both 
reducing the analytical and computational complexity of the representation and recovering the 
global ensemble dynamics. 

Inspired by the methods and models developed for crowd dynamics in I12[ 118) . in this 
paper we propose a new approach to the modeling of cell populations, based on measure-theoretic 
probabilistic tools, which permits a coherent representation of the system at both discrete and 
continuous scales as well as across intermediate levels. While the mathematical structures used 
here come to be very close to those firstly proposed in the above-mentioned works on crowds, the 
probabilistic point of view is instead original and finally makes the difference with respect to more 
traditional approaches, where the "discrete vs. continuous" dichotomy is referred to the spatial 
scale at which the system is modeled. In our approach, the underlying physical system (i.e., the 
cell population) is always regarded as a discrete (granular) one, thus all quantities can be directly 
referred to the scale of a single cell. However, we admit that the positions of the cells arc random 
variables. The key point is to describe the spatiotemporal evolution of the system by following the 
probability measure (law) attached to each of these variables rather than the variables themselves. 
This way we naturally pass from a Lagrangian to an Eulerian description, with the possibility 
to derive discrete or continuous models from the very same abstract framework by acting on the 
spatial structure of the probability distribution (atomic or absolutely continuous with respect 
to Lebesgue, respectively). Since the resulting models are such that the spatial distribution of 
cells affects the (nonlocal) dynamics of intercellular interactions, it turns out that the discrete or 
continuous character of the aforesaid probability can be related to the sensing ability of the cells 
themselves, thus providing quite naturally the concepts of multiscale dynamics and scale of the 
interactions. 

After the above general presentation, the contents of the paper can now be outlined more 
precisely. This Introduction is concluded by Section 11.11 which collects the main notations and 
theoretical notions used in the rest of the paper. In Section[5]the mathematical structures based on 
time-evolving probability distributions are introduced for modeling cell populations by a measure- 
theoretic approach. In particular, cell-cell interactions are discussed in some detail, including 
interactions mediated by a chemical possibly released by the cells themselves for signaling purposes. 
In Section [3] discrete and continuous models are addressed. A quick literature overview is offered, 
which shows that many recent microscopic/discrete and macroscopic/continuous models of cell 
colonies fit formally into our general framework. In addition, the latter can also help explain the 
links between the two levels of description, provided the convenient point of view of probability 
distributions is taken. Finally, the continuous version of our nonlocal model is studied in the limit 
of small sensing radius of the cells and local approximations, in the form of nonlinear degenerate 
diffusion equations, are derived and critically analyzed. In Section U the measure-theoretic tools 
are further exploited toward a hybrid approach, that introduces the above-mentioned concept of 
multiscale dynamics. In Section [5] a computational analysis of our method is performed by means 
of illustrative examples, which visualize the influence of the sensing radius and the sensing ability 
of the cells on the resulting colony dynamics. Finally, in Section [6] conclusions are drawn and 
possible developments of the method are outlined. The paper is also equipped with Appendix IA1 
which proposes a conceivable derivation of the probability measure equations out of underlying 
microscopic dynamics of individual cells. 

1.1 Notations and background 

In this section we quickly review some technical aspects, that will be extensively used in the sequel. 

Variables and nondimensionalization. We use the space x g M. d and the time t g [0, T max ], 
Tmax > being a certain final time (possibly +00), as independent variables. The space dimension 
d is 1, 2, or 3 for physical purposes. 

Throughout the paper, we will implicitly assume that these variables are dimensionless viz. 
referred to some representative values. Particularly, since a significant quantity is cell motility 
in free run conditions, the non-dimensionalization of space and time is related to the typical cell 
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Figure 1: Pictorial representation of the push forward via a mapping / in M. d . After the action of 
/, the measure of a set E is recovered by looking back at its pre-image under /. This way, the 
measure redistribution operated by / is conservative. 

speed, say V: if space is scaled with respect to a characteristic length L (e.g., the characteristic 
size of a bounded domain) then time is scaled with respect to L/V. As the "spontaneous" motion 
of a single cell is often regulated by diffusing chemicals, L and V might have a connection with 
the diffusion dynamics of the latter. For example, LV may be of the same order of magnitude of 
their diffusion coefficients. 

Given the reference values above, we will also tacitly assume that all other dependent variables 
and parameters are scaled accordingly. 



Functions. A function / of t, x is generically thought of as defined on [0, T max ] x R d . Instead 
of the classical f(t, x) we prefer the notation ft{x), which means that we are looking at / as 
parameterized by time (thus ft, for fixed t, is a function defined in M. d ). In particular, the subscript 
t will never stand for time derivative. 

The indicator function of a set A is 1a, thus 1a(x) = 1 if x € A, 1a(x) = if x $ A. 

The space of real-valued infinitely differentiable functions with compact support in R d is 
C c °°(R d ). 

Probability measures. We endow R d with the Borel tr-algebra £$(M. d ). We work with probabil- 
ity measures defined on B8(M. d ) and parameterized by time, that we denote by fit- More precisely, 
if <^P(R d ) is the space of probability measures on R d then fit is a mapping t £P(M. d ). 

The integral on M. d of a function / with respect to fit will be denoted by L d / dfit or, if it is 
necessary to specify the variable of integration, by J Rd f{x) dfi t (x). 

The push-forward of fit by a Borel-measurable function / : M. d — > R d is a conservative redis- 
tribution of fit in R d , which generates a new probability measure denoted by fifcfit and defined 
as 

(f#fi t )(E) := fi t (f-\E)), V£e^(R d ), 

see Fig. [TJ 



Random variables and stochastic processes. When dealing with random variables, we con- 
sider them as defined on an abstract probability space (f2, J^, P) and valued in the measurable 
space (R d , 33(M. d )). Here, is an abstract set (i.e., a set with no special structure), & is a a- 
algebra of subsets of J7, and P is a probability measure on & ' . For stochastic processes, we further 
assume that they are progressively measurable with respect to a filtration (=^t)te[o,T max l- 

The expected value of a random variable X : ft — > R d is denoted E[X]. Depending on the 
context, it can be computed by one of the following formulas: 

E[X] = J X{uj)dP{uj) = J xd{X#P){x), 
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X#P e ^(R d ) being the law of X. 



2 Modeling cell colonies by time-evolving probability mea- 
sures 

Following the ideas outlined in [21] , we model the position of a given cell in the physical space Mr 
as a random variable X\ valued in the measurable space (M. d , £%(W. d )). If the total number of cells 
in the colony is N then i = 1, . . . , N. However, instead of looking at the individual paths of the 
cells, we abstract their spatiotemporal evolution into that of the probability distribution attached 
to each X\. Therefore, the task is to devise a proper evolution model for such a distribution. 
A practical assumption for reducing the complexity of the problem is that cells within the same 
population are indistinguishable from one another. This is translated, in mathematical terms, as 
the X"s being identically distributed: fi t '■= X\^P is their common law at time t, independently 
of the index i labeling the cell. 

The key observation, leading to an equation for fi t , is that the knowledge of the probability 
distribution allows one to compute the (average) number of cells contained in a given set E C M. d 
at time t. Indeed, let Nf : 57 —> N be the random variable counting the cells in E € £$(M. d ) at 
time t, i.e., 

N 

N t E = j2^(xi)- 

i=l 

The aforesaid average number can be understood as the expected value of iV t s , which is 

N ^ N . 

mf] =J2J wK^dP^) = ]r y t E {x)dti t {x) = Nnt(E). (i) 

This simple formula matches the intuitive idea that the probability of finding a cell in a given 
set at a given time is the ratio between the average number (over several repeated experiments) 
of cells found in that set at that time and the total size of the cell population. ¥,[Nf] can be 
thought of as the mass of the set E at time t. In particular, from the relation (jXJ) we infer that the 
mapping E K[Nf] is, for fixed t, a measure proportional to fit, that we denote by m*. Since 
mt(M. d ) = Nfit(M ) = N is constant in time, i.e., the total mass of the system is conserved, it is 
coherent to assume that the measure m t satisfies the conservation law 

dm t „ , . 
-gf + V- (mt«t) = 0, 

where v t = v t (x) € R d is a suitable transport velocity. Substituting Eq. ([I]) gives an analogous 
evolution equation for the probability /Jt- 

^ + V • [tx t v t ) = 0. (2) 
Remark 1. Equation ^fy has to be understood in the weak sense of distributions, namely 

jfi;dfx t =f^-v t dfx t , V^eC c °°(M d ). (3) 

The interested reader is referred to [2Tj for a precise characterization of the concept of solution to 
such an equation. 
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2.1 Cell-cell interactions 



The dynamics of cell colonies are significantly determined by mutual interactions. This is because 
cells, as living matter, are not passively prone to Newtonian laws of inertia but are able to actively 
develop behavioral strategies depending on the presence of other neighboring cells and of chemical 
fields. To some extent, the modeling of cell dynamics can take advantage of certain ideas developed 
for other systems of interacting living subjects, such as human crowds and groups of animals, see 
e.g., jTOl [11] . Particularly, we borrow the idea that the active behavior of cells, as opposed to the 
passive one of inert matter, is translated into mathematical models by a direct modeling of the 
velocity, that cells are able to tune for themselves, rather than of the acceleration. 
The following general form of Vt, to be inserted in Eq. @, can be proposed: 



• T[ct] : R d — > K d is a taxisdike motion field generated by the spatial distribution of some 
chemical, whose concentration at time t is ct, that cells are sensitive to. We recall here, 
among other examples, the cases of chemotaxis, angiogenesis, vasculogenesis. T is a taxis 
operator, for instance for chemotaxis it is the spatial gradient, T[ct] oc Vet. This component 
of the velocity expresses indirect, or mediated, interactions among cells (particularly when 
the chemical is released by cells themselves for signaling purposes like in vasculogenesis). 

• K : M. d x R d — > M. d is the interaction kernel, expressing how a cell in x, hereafter called test 
cell, would interact with a neighboring cell in y, hereafter called field cell. Notice that only 
linearly additive binary interactions are taken into account in the present framework (see 
Section IQl for further discussions on this issue). This component of the velocity expresses 
direct mean field interactions among cells. 

A typical form of K is as follows: 



with k : R — > M. , i.e., the interaction between the test cell and a field cell depends on their 
relative position. When K has this structure, the interaction integral in Eq. (j4]) can be read 
as a convolution between k and the measure /it- 



In model direct interactions are nonlocal because the integral over M. d encompasses both 
near and far field cells. Indirect interactions may instead be either local or nonlocal, depending 
on the operator T ■ In both cases, they may have long-range effects due to the chemical dynamics. 
We refer the reader to Appendix [A] for motivating the form (Q} of cell velocity starting from 
microscopic cell dynamics written in terms of the random variables X\. 

A biologically reasonable assumption is that direct interactions extend over a finite range from 
the test cell. This is translated, in mathematical terms, by letting the function K(x, •) have, for 
all x £ R d , compact support, say 

U R {x) := suppi^z, •) C Br{x), 

where Br{x) C R d is the ball centered at x with radius < R < +00. Consequently, the interaction 
integral in Eq. Q can be rewritten equivalently as 




(4) 



where: 



K{x, y) = k(y- x), 



(5) 



N k(y- x) dnt{y) = N(k * ii t ){x). 




N 




Ur(x) 
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In the sequel, we will call R the sensing radius and Ur{x) the sensing neighborhood. Notice that 
the latter may either coincide or be a proper subset of Br(x). In the first case interactions are 
said to be isotropic (i.e., equally expressed all around the test cell), in the second case anisotropic 
(i.e., expressed only when external stimuli due to the presence of field cells come from specific 
directions). 

The chemical concentration can be assumed to obey classical reaction-diffusion dynamics, which 
implies the following prototypical evolution equation for Ct- 



dc t 
dt 



V • (£>Vct) + aN^ t - 704, 



(6) 



where D £ W dxd is the diffusion tensor (symmetric, positive definite), a > is the rate of 
autocrine production of the chemical by cells, and 7 > is the degradation rate of the chemical 
in the environment. 

Remark 2. Equation © is meant again in the sense of distributions: 
d 



dt 



ipc t dx= / V • (D T Vij)ct dx + aN ipdm-j t/jc t dx, Vip G C c °°(R d ) 



D T e R dxd being the transpose of D. 

By putting Eqs. ([2]), (j4|), ([6]) together, the following Eulerian model is finally obtained: 



dt 



dc t 



Ht T[ct](x) + N / K{x, y)dfx t (y) 



(7) 



— = V • (DVct) + aNiH - jet, 
k at 



A qualitative analysis of the initial- value problem associated with model (|Tp. concerning in 
particular the existence of measure-valued solutions t 1— > [J4 and the convergence to them of a 
properly conceived approximation scheme, is performed in the already cited paper [3T], however 
disregarding the coupling with the chemical field ct- Out of such an analysis, suitable forms of the 
interaction kernel K are discussed in both the isotropic and anisotropic cases. 

3 Discrete and continuous cell models 

The mathematical structure ((?]) can take advantage of the measure-theoretic formulation to con- 
stitute a paradigm for both discrete and continuous models of cell colonies. 



Discrete cell models. Discrete cell models are obtained from the assumption that the proba- 
bility /i t is atomic, i.e., 



1 N 



(8) 



h=l 



where X t := {xh (t)}h= 1 ^ s a time- varying set of N points in M. d . This amounts to saying that the 
random variable X\ is a discrete one, uniformly distributed in X t . In other words, X\ takes the 
values x\(t), . . . , XN{t) each with probability 1/N. In this case, the mass mt of a set E e ^(WL d ) 
is 



A' 



m t 



(E) - Nih(E) = ]T S Xh(t) (E) = card(i? n X t 



h=l 



i.e., the number of points xi(t), . . . , xw{t) contained in E (counting measure). Notice that, owing 
to the hypothesis of indistinguishability, it is not possible to associate deterministically the i-th 



G 



cell with its actual position in WL d at time t (it can be, with equal probability, any of the points in 
X t ). However, following the time evolution of the set X t is enough for an Eulerian description of 
the cell population as a whole. 

Plugging Eq. © into Eq. ^ gives 



jv 



h=l 



dx h (t) 
dt 



N 



J2 VV(ar fc (i)) • v t {x h (t)), W e C™(R d ) 



h=l 



whence — Vt{xh) in view of the arbitrariness of ip. Using Eq. Q for evaluating vt(xh) and 
considering the coupling with Eq. ^ ultimately yields 

N 



dx 



^=T[c t }(x h ) + J2K(x h , 



dt 

dc t 
dt 



Xl ) 



h = l,...,N 



1=1 



N 



(9) 



\7 ■(DVc t ) + aJ2 S x l (t)-lc t . 



1=1 



Model ^ is technically a hybrid instance of model ([7]) , as it is based on a discrete representation 
of the cell colony within a continuous description of the chemical field. 

Continuous cell models. Continuous cell models are obtained from the assumption that the 
probability p t is absolutely continuous with respect to the Lebesgue measure in Mr. Then there 
exists a nonnegative integrable function p t such that 



dfi t = p t dx, 



(10) 



or pt{E) = J E pt(x) dx for all E € 3§(R d ). In this case, the random variable X\ is a continuous one. 
Furthermore, also the mass m t is absolutely continuous with respect to the Lebesgue measure, its 
density being p t := N p t . The latter can be interpreted as the mass density of the cell population, 
whence the mass of a set E £ 3§(M. d ) is obtained as 

m t {E)=Nnt{E) = ( Pt {x)dx. 



By using (|10p in Eq. ([3]) , one recognizes the weak form of the conservation law 

^ + V • (p t v t ) = 0. 



With the velocity field (|4]) written for the measure (TT0|) : 

v t (x) = T[ct](x) + N K(x, y)p t (y) dy 



(11) 



and coupled with Eq. ©, this gives rise finally to the model 
dp t 



p t T[c t }(x)+N / K(x, y)p t {y)dy 



dt V ' 

= V • (DVc t ) + aNp t ~ jet, 

k at 



(12) 



namely a fully continuous instance of model ([7]) . 
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3.1 Quick overview of some pertinent literature 

At this point it is useful to compare the description of cell populations delivered by the measure 
equation 



Of 



7VV • 



tk / K(x, y)d^ t {y) 



= 



(13) 



(cf. Eqs. ((2J, Q decoupled from chemical dynamics, i.e., with T = 0) with the extensive lit- 
erature available on the topic. This can indeed contribute to a deeper understanding of our 
approach, which, as it should be clear from the previous sections and from Appendix [A] is con- 
ceptually different from taking directly either the microscopic or the macroscopic point of view. 
As discussed before, a formal comparison is possible with both microscopic/discrete and macro- 
scopic/continuous models. Therefore, for future reference, we explicitly write Eq. (|13[) with the 
discrete probability ([8]): 

x ]K(x h ,xi), h=l,...,N (14) 



df 



N 

E 



and with the absolutely continuous one (fT0|l: 



9p t 

at 



Pt / K(x, y)p t (y)dy 



= 0. 



(15) 



To begin with, it is particularly instructive to examine the paper |13) . in which the derivation 
of models at either scale is discussed through a progressive coarsening of the spatial resolution. 
Microscopic dynamics are assumed to be ruled by a Langevin's equation of the form 




N 

vf ) = J K m {t- t) Vi (t) dr + Fij + £00 , 



(16) 



(fi) 



(4) 



where m„ Vi are the mass and velocity, respectively, of the i-th cell (out of a population of N 
cells) and the various terms have the following meanings: 

(1) is the inertial force acting on the i-th cell; 

(2) is the (viscous) friction of the i-th cell with the outer environment, 7 > being the friction 
coefficient; 

(3) is the (viscous) friction of the i-th cell with the surrounding cells, in particular v? is the 
velocity of the i-th cell at the contact point with the j-th cell (and similarly for «V) and 
7ij > is the cell-cell friction coefficient; notice that this term accounts implicitly for the 
characteristic size of a cell; 

(4) is a memory term, with related memory kernel K m , expressing the tendency of a cell to keep 
the direction of motion (persistence) ; 

(5) is the resultant of pairwise actions exerted on the i-th cell by neighboring cells; 

(6) is a zero-mean time-uncorrelated random motion of the i-th cell. 

Conceivable simplifications of these general dynamics are obtained by neglecting inertia (1, friction- 
dominated motion), cell-cell friction (3), and persistence memory (4). Up to rescaling the noise as 
fi(t) = y/2j 2 Drji(t), where D > is a random motility coefficient, Eq. flTBl takes then the form 



dxi 
~dT 



1 



N 



F t3 +V2D m (t), 



(17) 
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where Xi is the position of the z-th cell (thus %^ — Vi). Equation ([!?]) resembles Eq. (fill but for 
the noise term, which is indeed possibly meaningless in an averaged representation such as the one 
delivered by probability measures (cf. Appendix \K§ . Actually, since Eq. (|17l) is meant to model 
cell-based microscopic dynamics, its counterpart in our approach is more properly Eq. (|26[) for the 
random variables X\, with which the analogy is evident if one sets T = 0. 

Similar ideas are used also in |15j , where an equation of the form (|17[) , though without stochastic 
noise, is derived under the assumption of friction-dominated dynamics. In this case, the formal 
comparison is with our discrete equation (1141) . which however describes actually the evolution of 
the support X t of the probability measure ©. 

In order to pass to a macroscopic model, in |13| an intermediate stochastic description by a 
cellular automaton is adopted, whose evolution rules are inspired by the microscopic model (|17[) . 
In more detail, the space is partitioned in a number of sites, in each of which cell division and 
apoptosis at respective rates A, C, > are taken into account, as well as hopping between adjacent 
sites at rate <f> > 0. By further coarsening the spatial resolution, the following equation for the 
cell density p t — pt(x) is finally obtained: 

dpt a ^ + 1 



-Ap t + [A(l - Pt ) - (]p t + vTkptVt, 



dt 2d 

where d is, like in our setting, the geometric dimension of the spatial domain and £ is the charac- 
teristic size of a cell. The former is a stochastic partial differential equation, which, in the absence 
of cell duplication and death (A = ( = 0), gives rise to a linear deterministic local model for cell 
motion: 

dpt P<P. 

to be compared with the nonlinear local models that we obtain from Eq. (1151) in the limit of small 
sensing radius (cf. Section [32]). 

Other authors [21 OH HZ] take directly a macroscopic point of view, modeling cell dynamics on 
a phenomenological basis by a conservation law of the form 

^ + ^-(p t F)=DAp t , (18) 

where D > is the diffusion coefficient showing up at the macroscopic scale as a consequence of 
random microscopic cell dynamics and F is the total force acting on a cell due to the presence of 
other surrounding cells. The coefficient cj) > is a kind of viscosity of the cellular matter seen as 
a continuum and R > is the sensing radius of the cells. Specifically, F is written as the sum 
of individual contributions within a sensing neighborhood extending up to a distance R from the 
cell: 

F = a f n{\y-x\)^-g{p t {y))dy, (19) 

J \y-x\ 

B R {x) 

see also [S]. Notice that ^F is the advection velocity of the cells, to be compared with our 
Eq. (fTT]) (with T= 0). In particular, if we assume that our interaction kernel is of type ([5]) with 
k compactly supported in the ball Br(0) then setting k(z) = afl(z)z/ \z\ and choosing g{z) = z 
makes Eqs. (fTS"]) and (flBl) formally consistent with one another (but for the linear diffusion term 
at the right-hand side of Eq. (|18[) , that however we can easily include in our equations as shown 
in Appendix \^ . Some comments about the relationships between our model (|15p and model 
([15 ]) - ([19" ]) are in order. 

• Interactions described by Eq. (|19p are, in general, nonlinearly additive whenever the function 
g is nonlinear. In contrast, our model features linearly additive interactions only. Concerning 
this, it is worth noticing that nonlinearly additive interactions are technically possible only 
when the density of cells can be defined. They cannot be straightforwardly extended to the 
case of more general mass measures (e.g., Dirac measures) because the action of a function 
g on an abstract measure is typically not defined. 
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• Model ([T5]) - ([T5t relies on the concept of sensing radius of the cells, which is however not rig- 
orously justifiable at the macroscopic scale because it pertains intrinsically to a microscopic 
length scale. Instead, with our approach it is fully acceptable also in a continuous context, 
since "continuous" refers here to the space structure of the probability distribution of cells, 
the underlying physical system remaining discrete and microscopic. 

Finally, we want to mention the chemotaxis model proposed in |14j , which describes cells as a 
continuum with density pt subject to interactions mediated by a chemical factor ct- 

^ +V • (xPtVrCt) = D p A Pt 
dct 

— = D c Ac t + apt - 7Q. 
at 

Here, V r denotes a nonlocal gradient sensed at a distance r > from the cell: 

9Bi(0) 

These equations can be duly compared with our continuous model (|12p by setting K = Q (and 
including the random motion of cells as explained in Appendix |A"]) . In particular, we notice that the 
nonlocal gradient may be a valuable hint for modeling the taxis operator T, considering that, for 
the same reasons discussed above, our framework enables us to handle rigorously the non-locality 
introduced by the parameter r as a microscopic feature of individual cells. 



3.2 Nonlocal vs. local continuous models 

Unlike those presented in Section [3.11 several other models of population dynamics are formulated 
assuming local (namely, pointwise) interactions, usually by means of nonlinear advection-diffusion 
equations at a continuous level [H [SJ [53]. In this section we investigate which kind of local 
continuous cell models can be obtained from the nonlocal equation (IT51) (cf. also [22]) under the 
assumption of small sensing radius. In particular, we discuss suitable analytical forms of such 
models in view of the locally dominant cell behavior (attraction or repulsion). For the sake of 
simplicity, we confine ourselves to the one-dimensional setting (d = 1). We consider interactions 
depending on the relative position of the cells (cf. Eq. ([5])), addressing separately the isotropic 
and anisotropic case. 

3.2.1 Isotropic interactions 

The sensing neighborhood of the point x is Ur(x) = (x — R, x + R), hence the velocity takes the 
form 

x+R 

v t (x)=N J k{y-x)p t {y)dy. (20) 

x-R. 

Notice that, in the interaction integral, it results \y — x\ < R. If R is small and we assume that k, 
p t are sufficiently smooth, we can expand the first one about zero and the second one about x as 

k(z) = fc(0) + fc'(0)* + ifc"(0)z 2 + ifc"'(0)z 3 + ^k IV {Q)z A + o(z 4 ) 

p t{y ) = p t{x ) + ^( X )(y -x) + \^{x){y - xf + ~0Or)(l/ - *) 3 (21) 

1 Pt {x)(y-xf + o{{y-xf), 



24 dx A 
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whence 



- ■>■),,,[!)) = A-(<l)/->,(.r) ! ( k{0)^(x) + k'(0)p t (x) j (,, - ,") 



+ 



^*(O)0(x) + fc'(O)^ (*) + ifc"(0)&(*)) (y - .x) 



+ ifc"'(0)p t (a:))( 2/ - a; ) 3 



+ i^, + >)g w+ ^4(,) 

+ ^"'(0)^ 0*0 + ^»^)) (y - *) 4 + o((y 



Now, since 




(y - x) n dy = 

" r I — Ji even, 

x H y n + 1 

we obtain from Eq. (f!?0"|). neglecting the term o((y — x) 4 ), 

= C p t (x) + C^{x) + C 2 ^(x) + C 3 ^(x) + C^(x), 
where the coefficients Cq , . . . , C4 are 

Co = (2k(0) + h"(0)R 2 + ±k IV (0)R 4 ^j NR 
Cl = Qfc'(O) + ^U"'(0)i? 2 ^ NR 3 

C 2 = Qfc(O) + ^fc"(0)i? 2 ) NR 3 

C 3 = ^-k'(0)NR 5 
15 

Ci = ±-k(0)NR 5 . 
60 

We can further specialize them by introducing the following hypotheses on k: 
(HI) fc(0) = 0, for excluding self-interactions 
(H2) k(—z) = —k(z), radial isotropic interactions 

which imply that even-order derivatives of k vanish in zero. Consequently: 



15 

?5 



Co 


= c 2 























hence Eq. (1151) is finally approximated as 



dp t 
dt 



n 9 ( d Pt \ , n d( d 3 p t 



OX \ OX I OX \ OX 
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Figure 2: Examples of interaction kernels k (cf. Eq. (0) for repulsion (left) and attraction (right), 
with sensing radius R. In the first case, field cells located in a left neighborhood of the origin 
make the test cell in z = move preferentially to the right, whereas those located in a right 
neighborhood make it move preferentially to the left. This is repulsion from populated areas. In 
the second case, the behavior of the test cell is the opposite, implying attraction toward populated 
areas. 



By introducing the rescaled space x* := R 1 x and time t* 
p, C i-i RH »(Rx*) =: pp(x*), we restate Eq. (231) as 



3fr 
dt* 



Pt* 



dp* t * 
dx* 







Ci\R 2 t and redefining pt(x) 

d 3 p* t . 



R 2 \d\ dx* 



dx* 



In the approximation C\ = Ifc' (0)NR 3 , which holds for small R, it results sign(Ci) = sign(fc'(0)) 



and ~~ mjbi i 



- sign (fe'(O)), thus we ultimately get 

dp*r 
dt* 



sign(fc'(0)) 



\—( 




Id/ 




dx* \ 


Pt ' dx* ) 


+ lOdx 7 V 


Pt dx* 3 )\ 



This equation combines second and fourth-order nonlinear degenerate diffusion. Its character 
is strongly determined by the sign of &'(0), in particular we can distinguish the following cases: 

• If fc'(O) < then sign (fc'(O)) = — 1, which implies that second-order diffusion is the classical 
forward one, whereas fourth-order diffusion gives the equation a backward character to the 
leading order. Since backward diffusion equations are ill-posed, especially in terms of con- 
tinuous dependence of their solutions on the data, when fc'(O) < it is advisable to stop the 
expansions ((2Tj) at the second order, so as to obtain a well-behaved diffusive local model of 
the form 



dt* 



d 

dx* 



Pf 



dp;, 

dx* 



Notice that this is consistent with the fact that fc'(O) < describes a locally repulsive 
(disaggregating) behavior of the cells (see Fig. (5J left). 

If fc'(O) > then the picture is the opposite to the previous one. Now sign(fc'(0)) = 1, thus 
second-order diffusion is backward and expansions (|21[) up to the fourth-order are necessary 
in order to obtain a formally well-behaved model at the leading order: 



dp;. 

dt* 



d 

dx* 



Pi- 



dp;, 

dx* 



i a 



d 3 p; 

dx*' 



Again, this is consistent with the fact that fc'(O) > describes a locally attractive (aggre- 
gating) behavior of the cells (see Fig. [2j right). 
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3.2.2 Anisotropic interactions 

In one space dimension, the prototypes of anisotropic sensing neighborhood are the half intervals 
Ur{x) — (x, x + R) and Ur(x) = (x — R, x). In this case, the coupling with the chemical dynamics 
can be exploited for determining the local orientation of Ur. For instance, the latter can be chosen 
in the direction of the taxis motion field in order to reproduce cell polarization. Hence the above 
alternatives apply for T[ct](x) > and T[ct](x) < 0, respectively. If T[c f ] changes sign then Ur 
will not be the same for all points x. 

For illustrative purposes, in this discussion we fix Ur(x) = (x, x + R)- After performing the 
same expansions as in Eq. (I21[) . we observe that now 



x+R 



(y - x) n dy 



R n+1 

n+l' 



thus the new coefficients Co, . . . , C4 to be used in Eq. (|22l) are 

Co 



C 2 
C 3 

c 4 



fc(0) + ^'(0) J R+ifc"(0) J R 2 + ^fc"'(0) J R 3 + I ^^(0) J R 4 j XB 



~fc(0) + ~k'(0)R + h"(0)R 2 + ^fc"'(0)i? 3 ) A If 
h(0) + h'(0)R+^k"(0)R^NR 3 

k m + ^o k ' m ) nr4 



1 
120 



k(0)NR 5 



Imposing that k satisfies hypotheses (HI), (H2) yields again that, in the above general expressions, 
fc(0) and all even-order derivatives of k computed in zero vanish. Nevertheless, this time all 
coefficients of the even-order derivatives of pt in Eq. (f2"2"j) do not vanish because of anisotropy: 



Cv = \(k'{Q) + ^k'"{Q)R 2 ) Mf 
Ci = Uk'(0) + ^k"'(0)R 2 ) xir 



Co = 



h>(0)NR* 



1 



C 3 = ^k\0)NR 5 
C 4 = 0. 



Correspondingly, the local equation obtained for the density is 



d Pt , r 9 . 2 , _ 



C — (- +c A (- d2 pt 

dx \ dx J dx V dx 2 



dx V dx 3 



which, after rescaling space, time, and density as in Section ^. 2.11 reads 



dp: 



d 



g£ + -sign(fc'(0)) — (p?) = -sign(fc'(0)) 



d 

dx* [ Pt * dx* 



dp;, 



3 d 



d 2 p* t , 



+ 8dx* l Pt * dx* 2 



1 d 

10 dx* 



Pt* 



d 3 Pt* 
dx* 3 
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Notice the appearance of the odd-order terms d x *(pt*) at the left-hand side and d x * (p|»9|«/5|») 
at the right-hand side, which are symptomatic of the anisotropy of the interactions. Once again, 
the sign of the highest-order derivatives, and consequently the appropriate order of the expansions 
(f2"Tj) , depends on the sign of k' (0) , namely on the local character of the interactions. Considerations 
analogous to those developed in Section ^. 2.1l can therefore be proposed: 

• If k'(Q) < (repulsive or disaggregating interactions) then a second-order expansion is 
convenient to obtain a formally well-behaved local model: 

a& 3 9 2 = d f d%. \ 

dt* 2 dx* [Pt * ' dx* \ Pt " dx* J ' 

• If A;'(0) > (attractive or aggregating interactions) then a fourth-order expansion is necessary 
in order to compensate for the resulting backward second-order diffusion: 

d%> , 3^,.^ = _JL (~ dpy\ _ 3 _d_ f _ i_ _d_ f d 3 p* t , \ 

dt* 2 dx* [Pt ' ' dx* \ Pt * dx* ) 8 dx* \ Pt ' dx* 2 J 10 dx* \ Pt * dx* 3 J ' 

4 Hybrid approach for multiscale dynamics 

The measure-theoretic formulation of model Q can be further exploited to come to a hybrid 
approach, in which discrete and continuous representations of the same cell population coexist 
and complement each other. Technically, this amounts to expressing the probability pt as a 
combination of an atomic and an absolutely continuous part: 

1 N 

p t = (l-u)-J2^ h (t)+up t C d , (24) 
h=l 

where C d denotes the Lebesgue measure in R d and u G [0, 1] is a parameter. The purely discrete 
and continuous cases, discussed in Section |3l can be recovered for u = and u = 1, respectively; 
in addition, a full range of intermediate representations is now accessible for < u < 1. 

It is worth stressing that when u ^ 0, 1 neither the sole set of points X t = {xh{t)}h=i nor 
the sole density pt provides the actual space distribution of the cells. According to Eq. (|24p . the 
complete picture results instead from a it-weighted joint contribution of the twos. In fact, the 
number of cells contained in a set E € 3§(R d ) at time t is 

m t (E) = (1 - u) card(£ n X t ) + uN J p t {x) dx, 

E 

thus for it ^ there might be a nonzero probability of finding a cell in a set E which does not 
contain any of the points of X t , whereas for it^l there might be a nonzero probability of finding 
a cell in a set E where the density p t is zero. 

With the measure (|24p the velocity field Q takes the form 

N 

v t (x) = T[ct](x) + (1 - u) K(x, x h (t)) +uN K(x, y)p t (y) dy, (25) 

which shows that at every point x £ IR d direct interactions depend simultaneously on the discrete 
and continuous spatial distribution of the cells. 

From the modeling viewpoint, this hybrid approach can be given various, possibly comple- 
mentary, interpretations. In the following we critically discuss some of them, having in mind the 
specific application to cell populations understood as living complex systems. As a general fact, we 
notice that such an approach is different from other hybrid approaches proposed in the literature, 
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Table 1: Parameters of the model used for the numerical tests of Section [5j 



Parameter 


Symbol 


Value 






Fig. El 


Fig. a 


Number of cells 


N 


2 


5 


Sensing radius 


R 


0.3 


0.6 


Cell sensitivity/scale of interactions 


u 


0, 0.1, 0.9, 1 


Final time 


T 

- 1 - max 


20 



in which the discrete and continuous representations refer to distinct parts of the system [5J 119) in 
a spirit similar to model ([§]). Here, it is the spatial distribution of one single cell population that 
is represented at an intermediate level in probabilistic terms. 

A first way of understanding the hybrid representation (j24|) is as a multiscale coupling between 
the "microscopic" scale of the set X t and the "macroscopic" scale of the density pt, with u playing 
the role of an interpolation parameter. This interpretation is brought from the multiscale model of 
crowd dynamics proposed in |12) , also based on a measure-theoretic approach to the conservation 
law ©. 

A different point of view, that we promote here, is instead not directly concerned with the above 
space scale dichotomy. As already stressed, the modeling structures presented in Section [2] refer 
to an intrinsically microscopic observation of the cell colony, with cell trajectories however blurred 
into the probability distribution p, t . Hence, the hybrid approach concerns the representation of the 
space distribution of cells as an abstract mathematical object, the look at the underlying particle 
system remaining an individual-based one. In particular, according to Eq. (|25D . a more discrete 
or continuous probability distribution influences the outcome of direct interactions among cells. 
In this sense, it is possibly appropriate to speak of multiscale dynamics based on the scale of the 
interactions. The parameter u can then be related to a sensing ability of the cells. A value of u 
close to (discrete representation) may characterize a population in which cells sense precisely 
other surrounding cells, whereas a value of u close to 1 (continuous representation) denotes a 
colony of cells that generically feel each other as a whole. This concept will be further discussed 
and clarified in the next Section [5l 

5 Numerical tests 

In this section, we use numerical simulations for exemplifying various dynamical scenarios gener- 
ated by the presented model according to different representations of the cell distribution. In order 
to focus genuinely on the effects of multiscale dynamics, we disregard the coupling with diffusible 
chemicals, namely we rest on Eq. (fT5|l with the measure pt given by Eq. (j2"4")l . 

Two-dimensional (d = 2) simulations are produced by means of a modification of the numerical 
scheme developed in [T2] (cf. also [TH1 [H] for error and convergence analysis) , specifically designed 
for the measure-theoretic framework with hybrid pt- The simulation setting is as follows (cf. 
Table[T]for a summary of the parameters). The spatial domain is il = [0, 4] x [0, 4], with an initial 
distribution of N — 25 cells compactly supported in the sub-domain £Iq = [1, 3] x [1, 3]. In more 
detail, the atoms and the density of the measure at time t = are chosen, respectively, as: 

x h{0) = I — - — , — - — I , hi, h 2 = 1, ■ • ■ , 5 
p (x) = - l Qo (x), 

hence the discrete and continuous parts of the distribution are initially uniform in fig. We model 
aggregation dynamics using an interaction kernel K of the form ([5]), with k(z) = zlg^g) (z), 
resulting in 

K(x, y) = (y - x)l BR{x) (y). 
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Cell-cell interactions are therefore isotropic. We consider two possible values of the sensing radius, 
R = 0.3 and R = 0.6, respectively. 

Figure |31 refers to the case R = 0.3 for different cell sensing abilities, u = 0, 0.1, 0.9, 1. In 
all cases, only partial aggregation of the colony occurs, indeed the final configuration consists 
in multiple isolated clusters. However, an important difference is immediately seen between the 
two limit cases u — (discrete dynamics) and u — 1 (continuous dynamics). In the former, the 
extremely localized cell sensitivity makes the starting configuration a steady one: because the 
sensing radius is smaller than the initial intercellular distance, individual atoms do not aggregate. 
In the latter, the distributed cell sensitivity (which means that the test cell may feel, within its 
sensing neighborhood, the presence/influence of the surrounding field cells besides their exact 
position) triggers instead aggregation, leading to the formation of 9 main clusters interspersed 
with 4 secondary ones. In addition, shadows of additional 4, even less relevant, clusters appear in 
correspondence of the corners of the original distribution. 

Adding a small continuous contribution to the discrete setting (u = 0.1, so that 10% of the 
distribution is continuous and 90% is discrete) changes dramatically the dynamics with respect to 
the purely discrete case. Most of the mass now concentrates in 4 clusters, including all atoms of the 
hybrid measure fit', a secondary cluster is formed in the middle, however it does not contain atoms. 
This is due to the fact that purely discrete dynamics are "fragile" , because they strongly depend on 
the precise locations of the cells. An even small uncertainty in the initial cell distribution, implying 
a mild spatially continuous component of the measure fit, can in fact result in qualitatively very 
different results. 

Conversely, adding a small discrete contribution to the continuous setting (u = 0.9, so that 90% 
of the distribution is continuous and 10% is discrete) does not appreciably perturb the dynamics 
with respect to the purely continuous case. The same 9 main clusters appear, where all atoms 
tend to concentrate, again interspersed with 4 secondary clusters without atoms. Actually, in the 
proposed simulation, at the final time a single atom remains trapped in the upper-right corner 
of the original distribution, however it corresponds only to a probability of (1 — u)/N = 0.004 of 
finding a cell there (in terms of mass, 10% of the average total mass of the colony). 

From the above-described realizations, we can infer that purely continuous dynamics are much 
more stable to discrete perturbations than purely discrete ones are to continuous perturbations. 
This is confirmed also by simulations with different values of u, for instance u = 0.7 (not shown). 
Nevertheless, by further decreasing u the granularity of the distribution begins to modify the 
continuous dynamics, indeed the aforesaid 9 main clusters switch to 4 main plus 1 secondary 
cluster for u = 0.1. 

Figure [4] shows similar numerical experiments with R — 0.6. In this case, the final configuration 
of the cell distribution is a single cluster for all values of u, meaning that the sensing radius is 
sufficiently large for a full aggregation of the colony to always take place. However, the snapshots 
at intermediate times display different transient dynamics of aggregation, confirming again the 
fragility of the purely discrete ones as opposed to the robustness of the purely continuous ones. 

6 Discussion 

In this paper we have introduced a mathematical method, based on conservation laws and measure 
theory, which is able to describe the spatiotemporal dynamics of cell populations. The proposed 
approach is intended to cope with the difficulty of providing an Eulerian ensemble representation 
of cell colonies starting from a Lagrangian individual characterization of a microscopic test entity. 
More technically, the method consists in modeling cell positions in the physical space as random 
variables and then in looking at the evolution of their probability distribution (law) rather than 
at that of the random variables themselves. This way, the underlying physical system is always 
regarded as a microscopic granular one but its Eulerian representation can be either discrete or 
continuous according to the specific structure of the probability. 

In our opinion, this method offers several conceptual advantages over classical microscopic and 
macroscopic models. First, it enables one to transfer unambiguously microscopic concepts into 
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Figure 3: Multiscale dynamics of aggregation with sensing radius R = 0.3. In each plot, the atoms 
Xh (t) and the level curves of the density p t are shown. 
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Figure 4: Multiscale dynamics of aggregation with sensing radius R = 0.6. In each plot, the atoms 
Xh (t) and the level curves of the density p t are shown. 
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a spatially continuous description of the cell distribution. This includes, for instance, sensitivity 
to chemical fields and to other cells within proper neighborhoods. Secondly, in order to obtain a 
global representation not necessarily focused on individual entities, our approach does not force 
the assumption that a cell population is a continuous matter. This makes the method suitable 
to describe both large colonies and small aggregates, for which the limit of an infinite number 
of cells may not be acceptable. Finally, it permits to embed rigorously the intrinsic granularity 
(discreteness) of cell populations into a continuous distribution, giving rise to multiscale dynamics 
that can be linked to the sensing ability expressed by the cells. 

In short, the key idea of the method consists in a change of viewpoint: from the physical 
positions of the cells to the probability measure attached to them. It is worth noticing that, 
once the latter is available, all important information on the system, such as mass distribution 
and spatial patterns, is directly recoverable without going back to the detail of the aforesaid 
random variables. Therefore, the latter mainly serve as conceptual tools for deducing coherently 
a transport equation of the probability law. 

Illustrative simulations have shown nontrivial consequences of the multiscale dynamics on the 
predicted pattern formation, thus implying that the sensing ability of the cells and the related scale 
of their mutual interactions are crucial aspects that our framework can include. For this reason, we 
expect our method to give an interesting contribution in the development of multiscale approaches 
to the modeling of cell phenomena, particularly in those problems in which the interplay between 
discrete and continuous descriptions can be linked to different abilities and functions of the cells. 

A From microscopic dynamics to an Eulerian representa- 
tion of cell colonies 

Model ©-(Hi was derived mainly heuristically. This appendix aims at motivating it starting from 
the level of individual-based dynamics. 

The hallmark of the method we are going to apply is that the total number TV of cells, though 
arbitrarily small or large, is in any case finite. No limit for N tending to infinity will be taken, 
for this would prevent one from dealing at the same time, and within the same mathematical 
framework, with discrete, continuous, and hybrid representations. This is a significant difference 
between our point of view and other approaches, based on similar ideas, that have been used in 
the literature for bridging models at different scales in the limit N — » oo, see e.g., [SI 171 120]. 

A.l Cell-based microscopic dynamics 

In Section [2] we modeled the positions of the cells in the physical space M. d at time t as im- 
measurable random variables {Xj}^L l5 defined on the stochastic basis (fi, # ', (^t)t, P), with 
values in the measurable space (M. d , 3§(Mr)). We assume now, in more detail, that (X\ )te[o, r max ] 
is, for all i, a continuous-in-time process obeying some stochastic evolution equation. A prototype 
for chemical-induced taxis and binary interaction-based dynamics is 

dX* = \ T[c t ] (XI ) + J2 K ( X t . x t ) j dt + dWl , i = 1, . . . , N, (26) 

where a > is a constant (in the simplest case) and (Wj )te[o, T max ] are independent standard 
Wiener processes modeling stochastic fluctuations in cell motion (i.e., they summarize all aspects 
of cell behavior on which there is no direct phenomenological control) . 
A few remarks about Eq. (|26p are in order. 

• It describes the system from a microscopic Lagrangian point of view. In fact, cells are mod- 
eled individually and their paths t n- Xl(ui), ui S J7, corresponding to different realizations 
of the stochastic processes {Xl}f =l , are trackable in the physical space M. d . 
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• It models cell dynamics by controlling directly the velocity rather than the acceleration. 
This is consistent with the active behavior of cells, which, as mentioned in Section [2 are 
able to express strategies rather than being passively subject to inertia. 

• It assumes deterministic cell-cell interactions, in fact the latter are included in the determin- 
istic part of the right-hand side. This is reasonable for applications in which a standard (say, 
"rational") behavior of individual cells can be identified, over which very large deviations 
are not expected. 

A. 2 Ensemble representation by probability distributions 

We now pass to an Eulerian representation of the spatial distribution of the cells by means of 
probability measures. For this, let us pick a function ip G C^°(Mr) and compute it along the 
trajectories of system (|2"o) . From Ito's calculus it follows dif)(Xl) = V?/>(JQ) ■ dX\ + aAip(Xl) dt, 
whence 

t N t 

4>{X\) = ^{Xl) + J Vi>{Xl) ■ T[c t }(X l s ) ds + Y^J V^{X\) ■ K(X\, X{) ds 
o J' =1 o 

t t 



+ a J Aip(Xl) ds + V2<7 J Vip(Xl) ■ dW l s 
o o 

for all t € [0, T max ]. Taking the expectation of both sides and interchanging properly the order of 
integration at the right-hand side by Fubini's Theorem gives 

t 

il>(X*(u))dP(u) = J i>(Xi(u,))dP(u,) + J J Vilj(Xi(oj))-r[c t }(Xl(Lo))dP(cj)ds 
n n an 

N 1 



J=1 o n 



t 



a \ \ Aip(Xi(u))dP(u;)ds. 



o n 



Particularly, we notice that E[f * Vtp(Xl) ■ dWl\ — for all t € [0, T max ] due to the properties of 
Ito's stochastic integral. 

Let now n\ = X l t #P G 0>(R d ) be the law of X\ and fif = {X\, X\)#P G ^(R d x R d ) the 
joint law of X l tl X\. The previous formula can be rewritten as 

t 

ijj(x)dni(x) = J ^{x)d^{x) + J J V^(x)-T{c t }(x)dnl(x)ds 

R d 

N / * 

J J Vi>(x)- K(x, y) dtf {x, y)ds + a J J A^{x) d& (x) ds. (27) 



R d xR d R d 



In order to proceed toward an equation for the marginal distribution fj^, we disintegrate /i t J on 
R d x K d using the canonical projection n : R d x R d — >• R d , n(x, y) — x. Notice that fi\ — tt^A 1 
by definition and that, for all x G R d , the set 7r _1 ({a;}) C R d x R d can be identified with R d itself. 
Hence, owing to the disintegration theorem Q], for all i > and all j = 1, . . . , N there exists a 
family of probability measures {[ii ('\ x )}xe$t i <- ^(R d x parameterized by x, such that: 
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• iJ} l (-\x) is unique for /j,\-a.e. x G M d ; 

• supp/4'*(-|a0 C 7r _1 ({^}) - M d for /zj-a.e. x G M d ; 

• ${dxdy) = i4 {i (dy\x) <g> jtij(dx) for juf-a.e. (x, j/) G M d x M d . 

We notice that the second condition means that /ij "lives" actually on R d and can therefore 

be identified with an element of 3^{R d ). The third one is instead a restatement of the conditional 
probability theorem, up to understanding ^ l (-\x) as the distribution of X\ conditioned to X\. 
Particularly, it implies that integrals on K d x W l can be computed as iterated integrals on R d , 
thus we have: 

VtA(x) • K(x, y)d^(x, y)=J V4>(x) ■ ( / h'U. ,,) ,////''( ,,|.r) ) ,/;/j.ri. 

R d xR d R d 

Consequently, Eq. (|2"T)) becomes 

^(x)djuj(a;) = / ip(x)di4,(x) 



Rd 



W(x). (t[c 4 ](z)+^ f K(x, y) dfjpj i (y\x) j d/4(:r)<fe 



+ ct / / Aip{x) d^ s {x) ds 

R<1 

or, defining the velocity field 

:= T[c t ](x) + ]T j y) d^G/lz) (28) 

and denoting by (•, •) the duality pairing between distributions and test functions in C£°(M. d ), 

j t (f4, 1>) + (V ■ (nX), $ = cr(A^j, VtA G C c °°(K d ), a.e. * G (0, T max ] 

which is the functional form of the following advection-diffusion equation: 

V-(l4vi)=vAi4. (29) 



d(4 

— -r V .wtJ-»">f 

Equations (|28l) , (|29jl provide a deterministic Eulerian representation of the cell colony with cell 
trajectories t <-> X l t (uj) blurred into the spatiotemporal evolution of the probability distribution 
of cell positions in M. d . 

A. 3 Toward a self-contained modeling framework 

We discuss now a few manipulations that allow one to recover model ©-Q out of the general 
structure delivered by Eqs. (|28p. (|2"9"|) . and we critically analyze their validity and significance for 
modeling. 

We recall that in Section [2] we made the further assumption that cells are indistinguishable. 
In the present context this corresponds to the ansatz fij = ■ ■ ■ = /ij on the solutions to Eq. (|2U)) , 
which, for coherence, calls for the ansatz that also the conditioned law fi{^(-\x) does not depend 
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on the specific pair of indexes (i, j). Let fi t be the common law of the X"s and fi t (-\x) their 
pairwise joint law. Equation (l28l) takes the form 



while Eq. (|29l) specializes in Eq. © up to technically setting a = (namely, disregarding the 
effect of random walks of the cells) . 

With the further approximation fit(-\x) ~ /it for /x t -a.e. x £ M. d , Eq. (I3U1) becomes exactly the 
velocity field (£[]). We point out that the correct way of understanding such an approximation is 
in terms of the Wasserstein distance between /j, t (-\x) and fit in the space of probability measures 
with finite first moment. Indeed, it is this metric that is the most suited for model ©-(HI, cf. 
[21] , We also observe that the accuracy of this approximation may depend, among other things, 
on: 

• the assumed microscopic dynamics of the interactions, namely on the specific choice of K 
(for instance, if K = then the approximation is actually exact); 

• the total number of cells in the population. Heuristically, the larger N the lower the ex- 
pected statistical dependence between single pairs of cells, because the dynamics of a cell 
are presumably affected by many other neighboring cells. 

More accurate approximations, involving the assessment of the correlation between pairs of in- 
teracting cells, may be envisaged, adapting to the present context some arguments proposed for 
similar issues in different fields (cf. e.g., [TB] about vehicular traffic). 

From a different point of view, replacing // t (-|x) with /i t in Eq. (13T))) can be interpreted as a 
further modeling of the interactions among cells, that would be technically more difficult to im- 
plement in an individual-based framework. In more detail, it implies that the interaction strategy 
of the test cell does not depend on the past history: at each time, it reacts at field cells as if they 
had just been scattered randomly and independently in the surrounding space. Such a behavior 
can be justified in view of the indistinguishability of cells and of the non-cooperativeness of their 
dynamics (i.e., there is no jointly organized strategy among them). Detailing such an effect in 
Eq. (|26[) is perhaps not immediate and, after all, also needless, since the final interest is in the 
Eulerian framework ([2"51) -(12"§1). 

Finally, it is worth noticing that, in the case of a purely discrete representation of the cell 
population, the replacement of /i*(-|x) with /i t is fully consistent with Eq. ([25)1 . Indeed, if \it 
is given by Eq. (0) then the velocity (H]), evaluated at the h-th atom of \x tl coincides with the 
deterministic part of Eq. ((26|) . 

References 

[I] L. Ambrosio, N. Gigli, and G. Savare. Gradient flows in metric spaces and in the space of 
probability measures. Lectures in Mathematics ETH Zurich. Birkhauser Verlag, Basel, second 
edition, 2008. 

[2] N. J. Armstrong, K. J. Painter, and J. A. Sherratt. A continuum approach to modelling 
cell-cell adhesion. J. Theoret. Biol, 243(1):98-113, 2006. 

[3] N. J. Armstrong, K. J. Painter, and J. A. Sherratt. Adding adhesion to a chemical signaling 
model for somite formation. Bull. Math. Biol., 7I(l):I-24, 2009. 

[4] M. Bertsch, M. E. Gurtin, D. Hilhorst, and L. A. Peletier. On interacting populations that 
disperse to avoid crowding: The effect of a sedentary colony. J. Math. Biol., 19(1) : 1—12, 1984. 

[5] M. Bertsch, M. E. Gurtin, D. Hilhorst, and L. A. Peletier. On interacting populations that 
disperse to avoid crowding: preservation of segregation. J. Math. Biol, 23(1):1— 13, 1985. 




(30) 



22 



M. Bodnar and J. J. L. Velazquez. Derivation of macroscopic equations for individual cell- 
based models: a formal approach. Math. Methods Appl. Sci., 28(15):1757-1779, 2005. 

V. Capasso and D. Morale. Asymptotic behavior of a system of stochastic particles subject 
to nonlocal interactions. Stock. Anal. Appl, 27(3):574-603, 2009. 

V. Capasso and D. Morale. Stochastic modelling of tumour-induced angiogenesis. J. Math. 
Biol, 58(l-2):219-233, 2009. 

M. A. J. Chaplain, M. Lachowicz, Z. Szymahska, and D. Wrzosek. Mathematical modelling of 
cancer invasion: the importance of cell-cell adhesion and cell-matrix adhesion. Math. Models 
Methods Appl. Sci, 21(4):719-743, 2011. 

E. Cristiani, P. Frasca, and B. Piccoli. Effects of anisotropic interactions on the structure of 
animal groups. J. Math. Biol, 62(4):569-588, 2011. 

E. Cristiani, B. Piccoli, and A. Tosin. Modeling self-organization in pedestrians and ani- 
mal groups from macroscopic and microscopic viewpoints. In G. Naldi, L. Pareschi, and 
G. Toscani, editors, Mathematical Modeling of Collective Behavior in Socio- Economic and 
Life Sciences, Modeling and Simulation in Science, Engineering and Technology, pages 337- 
364. Birkhauser, Boston, 2010. 

E. Cristiani, B. Piccoli, and A. Tosin. Multiscale modeling of granular flows with application 
to crowd dynamics. Multiscale Model. Simul, 9(1):155-182, 2011. 

D. Drasdo. Coarse graining in simulated cell populations. Adv. Complex Syst., 8(2-3):319- 
363, 2005. 

T. Hillcn, K. Painter, and C. Schmeiser. Global existence for chemotaxis with finite sampling 
radius. Discrete Contin. Dyn. Syst. Ser. B, 7(1):125-144, 2007. 

P. Macklin, M. E. Edgerton, A. Thompson, and V. Cristini. Patient-calibrated agent-based 
modelling of ductal carcinoma in situ (DCIS) I: Model formulation and analysis. J. Theoret. 
Biol, 2011. (Accepted pending revision). 

P. Nelson. A kinetic model of vehicular traffic and its associated bimodal equilibrium solutions. 
Transport Theor. Stat, 24(l-3):383-409, 1995. 

K. J. Painter, N. J. Armstrong, and J. A. Sherratt. The impact of adhesion on cellular 
invasion processes in cancer and development. J. Theoret. Biol, 264(3):1057-1067, 2010. 

B. Piccoli and A. Tosin. Time-evolving measures and macroscopic modeling of pedestrian 
flow. Arch. Ration. Mech. Anal, 199(3):707-738, 2011. 

M. Scianna, L. Munaron, and L. Preziosi. A multiscale hybrid approach for vasculogenesis 
and related potential blocking therapies. Prog. Biophys. Mol. Biol, 2011. In press (doi: 
10 . 1016/j .pbiomolbio . 2011 . 01 . 004). 

A. Stevens. The derivation of chemotaxis equations as limit dynamics of moderately inter- 
acting stochastic many-particle systems. SIAM J. Appl Math., 61(1):183-212, 2000. 

A. Tosin and P. Frasca. Existence and approximation of probability measure solutions to 
models of collective behaviors. Netw. Heterog. Media, 2011. In press (arXiv: 1012.2395). 

S. M. Wise, J. S. Lowengrub, H. B. Frieboes, and V. Cristini. Three-dimensional multispecies 
nonlinear tumor growth - I: Model and numerical method. J. Theoret. Biol, 253(3):524-543, 
2008. 

T. P. Witclski. Segregation and mixing in degenerate diffusion in population dynamics. J. 
Math. Biol, 35(6):695-712, 1997. 



23 



